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We consider the operator product expansions for unitarity fermions. We compute 
the dynamic structure factor S(c[, ui) at large frequency and wavenumber away from 
the one-particle peak. The overall normalization of 5(q, u>) is determined by Tan's 
contact parameter, and the dependence on q and uj is obtained in closed analytic 
form. We also find energy deposited into the system by small, rapid variations of 
the inverse scattering length. 
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I. INTRODUCTION 

One of the most interesting systems currently under study is that of fermions with inter- 
actions fine-tuned to unitarity [1]. This system is interesting due to the absence of any scale 
parameter other than the density. The system is strongly interacting and thus presents a 
challenge for many-body methods. 

Recently, these systems have been probed with high-frequency or short-wavelength exter- 
nal probes. By bombarding the fermions with radio-frequency (RF) radiation, which causes 
a transition from one atomic state to another, one can probe the structure of the unitary 
Fermi gas. Previous investigations [2-4] have shown that the high frequency tail of the RF 
line shape is related to the contact parameter, first introduced by Shina Tan [5, 6]. 

In this paper, we address two problems. The first one concerns the dynamic structure 
factor S(q,u). This factor, in principle, can be measured by Bragg scatterings [7]. At 
the Fermi momentum and energy scales, the dynamic structure function depends on the 
complicated many-body physics of the unitary Fermi gas. For much larger q and u, the 
picture becomes simpler. There is a dominant coherent peak in the response function at 
oo = e q . Far away from this peak, the strength of the dynamic structure factor should be 
proportional to Tan's parameter, the same parameter that characterizes the RF tail. It is 
easy to understand why this is true: in order to absorb an external "virtual photon" carrying 
large energy and momentum and far away from the coherent peak, the photon has to hit a 
particle when it is near a second particle, so that it can give its momentum and energy to 
both. The absorption rate, therefore, depends on the probability of finding two particles at 
short distance from each other — which is characterized by Tan's contact parameter. 

The problem of the dynamic structure factor is more complicated than that of the RF 
response. In particular, care should be taken to not violate conservation laws. In this paper 
we use the operator product expansion (OPE) to facilitate the calculation. The operator 
product expansion, introduced by Wilson, is a standard method of quantum field theory. 
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The use of the OPE for cold atoms was pioneered by Braaten and Platter [8]. The end result 
for the dynamic structure factor is Eq. (26). 

The second problem considered in this paper is the calculation of the energy deposition 
into a unitary Fermi gas by small, rapid variations of the inverse scattering length. It can 
be solved using the same OPE methods. 

The usefulness of the operator product expansion can be illustrated as follows. Suppose 
we need to compute the following Green's function 

G AB (uj,q) = Jdtdxe^-^(A(t,x)B(0,O)) (1) 

for large u and q. Here "large" means energy and momentum much larger than the typical 
energy and momentum scales of the state with respect to which the average (...) is taken. 
For the ground state of a unitary Fermi gas, these scales are the Fermi energy and Fermi 
momentum. Let us also recall that one can associate a local operator O with a scaling 
dimension A<> In our counting scheme, the dimension of momentum is 1 and of energy is 
2 (the particle mass m is set to 1). Assuming that the operator product expansion exists, 
one can expand the product A(t,x)B(0 : 0) in terms of local operators, 

A(t, x)B(0, 0) = W A - A -- As /, O^Pj O,(0). (2) 

Here fa are functions of one variable |x| 2 /t, and A; are the dimensions of Oj. In contrast 
to the OPE in relativistically invariant theories, in nonrelativistic theories the OPE coeffi- 
cients are not constant, but are in general functions of this variable. This identity is to be 
interpreted as an operator identity; in particular, we can take the expectation value of both 
sides with respect to any state, including thermodynamic states. Taking the average and 
performing a Fourier transform, one finds 

Gab(u, q) = ^(s+^-L-a^^ Ojj) <°«>. 9= 111- ( 3 ) 

On the right hand side, the higher the dimension of 0«, the more rapidly its contribution 
decays in the limit of large momentum/energy (to be precise, the limit considered in this 
paper will be u — >■ oo, q — > oo, q 2 /u = fixed.) Thus, the leading behavior of the Green 
function is dominated by those few operators in the OPE with smallest scaling dimension. 

The expectation values of the operators Oj, are, in general, not computable theoretically 
because they depend on many-body physics. Thus, they should be considered as numbers 
parameterizing the many-body state. The OPE coefficients q, however, depend only on 
few-body physics (although, the number of bodies increases with increasing complexity of 
the operator Oj), and hence can be computed reliably, at least for simple operators A, B 
and Oj. In this way, the functional dependence of Gab on frequency and wavenumber can 
be expressed in terms of a few numbers which have to be determined experimentally or 
numerically. 
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We now discuss the operators of lowest dimension. The unitary fermions are described 
by the Lagrangian 

c = i^d t i> - ^£ + ^0* + 44(1) - c y 0. (4) 

We set the fermion mass to one. The indices 1, 2 refer to the two spin polarizations. If one 
integrates out 0, the saddle point for this field is at = Coip 2 i J i- This Lagrangian is there- 
fore equivalent to that with a four- Fermi interaction coipi^l^P^i- We will use dimensional 
regularization, where setting Cq 1 = corresponds to fine-tuning the interaction to infinite 
scattering length. The propagator of the field is completely determined by its self-energy, 
and is equal to 

£>(o;,q) = — =^=. (5) 
vV/4 -u-iO 

Let us discuss the operators which may have nonzero expectation value in a unitary 
Fermi gas. For our applications, we will need to consider only operators which do not carry 
particle numbers. Moreover, we assume the ground state to be isotropic, so we need to look 
only at operators with vanishing orbital angular momentum. The following two operators 
have the lowest dimensions: rii = ^fyi i} = 1 , 2) and 0*0. The operators rii are the particle 
number densities and have dimension 3. The operator 0*0 has dimension equal to 4. This 
can be seen by computing the dimension of from Eq. (5). An alternative way to find the 
dimension of is to use the operator-state correspondence, according to which the dimension 
of is the ground state energy of a system of one spin-up and one spin-down particle in a 
harmonic trap with unit oscillator frequency. This ground state energy is 2. 

As we shall see, the expectation value (0*0) can be identified with Tan's contact param- 
eter. Together, these two operators dominate the high-momentum behavior of correlation 
functions, including the dynamic structure factor. 



II. SINGLE-PARTICLE GREEN'S FUNCTION 



To establish the relationship between (0*0) and Tan's contact parameter, let us apply 
the method of the OPE to the one-particle Green's function 

iG(u,q) = J dtdxe iujt - i(iX <7Vi(f,x)^!(0,0)), (6) 

where T denotes time ordering. 1 At the end we will be interested only in the case of 
t < 0, where the Green's function is — (■0 t (O)'0(t, x)), but the Feynman diagrams are readily 
available for the time-ordered Green's function. 

According to the previous discussion, we can write 

G(u,q) =C n (a;,q)(^2) + C^(a;,q)(0*0) + --- (7) 



Y. Nishida [9] has essentially arrived at a similar derivation. 
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(We shall see why there is no (VqV'i) term below.) 

To compute C n , we use the Feynman diagram in Fig. 1. 




FIG. 1: The diagram that gives C n in Eq. (7). The large momentum and energy flow along lines 
carrying label "q." 

Physically, this diagram describes the interaction of a particle carrying large momentum 
and frequency (u, q) with particles already in the medium. The energy and momentum 
of the particles in the medium can be neglected compared to (cu,q). The "hard" particle 
interacts with particles of type 2 in the medium with the result, 



One can see immediately that C n (u,q) has singularities only in the lower-half plane, and 
thus does not contribute to the Green's function for t < 0. This can already be seen from 
the Feynman diagram in Fig. 1: one can go directly from the initial point to the final point, 
following the direction of the propagators (fermions and scalars). Recall that the propagators 
are all retarded. 

Consider now the OPE coefficient of C^*^. It can be found by computing the fermion 
propagator in the background field of 0. One diagram that contributes to C^*^ is as in 
Fig. 2. The dashed lines carry zero energy and momentum, and go to the "condensate" 



(0*0). Fig. 2 is not the sole diagram contributing to C<^; there is an infinite number of 
diagrams contributing to this coefficient, like those depicted in Fig. 3 and obvious subsequent 
iterations. It can be seen that, in all diagrams of the type of Fig. 3, one can go from one end 
of the diagram to the other by following the arrows, which is not true in the case of Fig. 2. 



C n (u,q) 



(u - q 2 /2 + iO) 2 v /q 2 /A-u - iO ' 



(8) 




FIG. 2: A contribution of (0*0) to the fermion propagator 
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(t,x) 
• — 



(0,0) 



FIG. 3: The other diagrams contributing to the fermion propagator 

Thus, if our ultimate goal is to compute the Green's function for t < 0, it is sufficient to 
just evaluate the diagram in Fig. 2. 

From the Feynman rules, the contribution of Fig. 2 to is 

1 



(w-e q + i0) 2 (-w-e q + i0)' 



(9) 



where e q = q 2 /2. 

With this information we can find the leading nontrivial contribution to the Green's 
function for t < 0: 



/ 



du 



dxe^ x (v4(0,0)<Mt,x)> = -i l^e-**Cp4(u,q)(<F<l>) 

1 exp (i^-t) (0*0) + 



In the limit t — > — 0, we find the tail of the distribution function to be 

(0*0) 



n n = 4 
q 4 



t < 0. (10) 



(11) 



which allows us to establish (0*0) = C, where C is Tan's contact parameter [5, 6]. 

One can also show that the diagrams in Fig. 3 do not contribute to the imaginary part 
of G(u,q) for uj < 0. The latter receives only a contribution from Fig. 2, which is 



u < : ImG(w, q) 



-7iS(co + e c 



(12) 



The peak at u = — e q is already discussed in Refs. [11, 12]. At this level, we are not able to 
resolve the structure of the peak. 

III. RF SPECTROSCOPY 



Consider now a system where, in addition to the 'up' and 'down' fermions included in 
the earlier Lagrangian, a third species of fermion is added which does not interact with the 
other two fermions. Suppose we turn on a photon field which converts atoms of type one 
into atoms of type three. The absorption rate is proportional to the imaginary part of the 
Green's function of 



O 



13 



(13) 
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We thus need to compute the OPE expansion of 

n(w,q) = -i I dtdxe™-*** (TO 13 (*,x)O[ 3 (0,0)> = C n (m> + 0^(0*0) + 



The coefficient C„ is obtained from the diagram of Fig. 4, 

Cr?. — 



(14) 



(15) 



uj — e q + iO 

The external wavy lines in this diagram represent the operator O13. The interpretation 
of this formula is obvious. The dominant part of the response function is the same as a 
collection of noninteracting particles with zero momentum, and has a peak at u — e q . 

3 

wxaaa/ Jwvw 

FIG. 4: Diagram contributing to C n 

Let us now compute the diagram which determines the coefficient of 0*0 in the OPE 
(Fig. 5). We find 

'dp dp 1 1 1 



n w,q) 



= ... + 



(27r) 4 (p - e P + ?0) 2 p + oj - e p+q p + e P - «0 

+ (a; ->■ -00, q ->■ -q) H 

dp 1 1 



(16) 



» + (w ->■ -u, q ->■ -q) + 



(2tt) 3 4e2 w - e p - e p+q + iO 
This integral is divergent in the infrared. We can isolate the divergent part of the integral 



as 



(0*0) 



dp 1 f dp 1 



»• (17) 



w - e q + iO J (27r) 3 4e£ J (27r) 3 4e p \w - e p - e p+q + iO co - e q + iO, 

It is clear that the first (divergent) part is actually included in the first term in the OPE 
(15), i.e., it represents the contribution of the condensate (0*0) to the number density n. 
Subtracting this contribution, the second part gives C^, 

dp 1 / 1 1 



Crh* r , 



(27r) 3 4e^ \uj — e p — e p+q + iO u — e q + iO 
We now restrict ourselves to the imaginary part at q = 0. We have 



+ (u ->■ ->• -q). (18) 



A „) = -I to „ (u ,0) = ^, 



(19) 



which coincides with the result obtained previously by other authors [3, 4]. 
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FIG. 5: Diagram contributing to C^*^. The external dashed lines, which carry zero energy and 
momentum, go to the "condensate" (0*0) which is not depicted for simplicity. 

IV. DYNAMIC STRUCTURE FACTOR 



Having tested the OPE method on the simple examples above, we will now discuss the 
dynamic structure factor. The dynamic structure factor can be defined as the imaginary 
part of the Green function of the density, 

G nn (u, q) = -ij dt dq e l "- x (Tn(t, x)n(0, 0)) = C n (n) + 0^(0*0) + • • • (20) 

Namely, 



%w) = --ImG n „(w,q) 1 u > 0. (21) 

7T 



The coefficient C n is 



C n = G(u, q) + G(-w, -q) = 2 2£ 2 \. n , (22) 

and corresponds to the one-particle peak ~ 5(u> — e q ) in the dynamic structure factor. This 
peak dominates all sum rules. However, we are interested in the structure factor far away 
from this peak, hence we need to compute the OPE coefficient of the next operator, 0*0. 

The diagrams contributing to C^*^ are sketched schematically in Fig. 6. 

The contribution from the first diagram is exactly twice the diagram in the RF case, 

C - = 2 / Wf 4 CVw>0 - ^T») (23) 
The second diagram gives 

c m = 2 f_*P L_ ( 1 m 

7 (2tt) 3 4e p e p+q Vw-e p -e p+q + iO u> + e p + e p+q - iOJ 

The third diagram is the product of three pieces, each of which can be computed separately. 
At the end one finds 



^(3) 



-2 



dp 1 1 



(27r) 3 2e p u - e p - e p+q + iO 



—An 

+ (u ->■ -u;,q ->■ -q). (25) 



|e q - u — iO 






FIG. 6: Schematic diagrams contributing to G$$. Each diagram represents a number of Feynman 
graphs where the photon lines are attached to different fermion lines 

Evaluating the imaginary part of the correlation function for positive u, we find that it is 
zero when u < q 2 /4. This should be the case since u = q 2 /4 is the threshold for knocking 
out two particles from the medium (recall that all the energy scales associated with the 
medium, like the chemical potential, are negligible). For u > g 2 /4, the dynamic structure 
factor is 



S(q,u) 



1 



ImCU 



1 " <?74 + 1 



7T 



2tt 2 (lo - q 2 /2) 2 2ir 2 ujq 
1 



In 



oj + q^Ju - q 2 /4 
\io-q 2 /2\ 



7T" 



Yvw-g 2 /4 



ln 2 " + ^-^/4 _ 7r2%2/2 _ u;) , (26) 



\uj-q 2 /2\ 



In Fig. 7 we plot the function u 3 ^ 2 S(ui, q) / as a function of the "Bjorken x" variable, 
x = q 2 /2u. This can be thought of as a plot of S at fixed u as a function of q 2 . In Fig. 8, 




0.5 1.0 1.5 2.0 

FIG. 7: A plot of w 3 / 2 S(q, u) /{</>*</>) as a function of x = q 2 /2u. 
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we plot the same quantity function of x 
as a function of u at fixed q. 



2ui/q 2 . This plot is basically a plot of S 



20 - 



10 




0.5 



1.0 



1.5 



2.0 



FIG. 8: A plot of g 3 / 2 5(q,o;)/ (0*0) as a function of x" 1 = 2oj/q 2 



As we see from the plots, the dynamic structure factor is positive for all u and q, as it 
should be. It vanishes in the zero momentum, finite frequency limit q = 0, as it should 
due to the conservation of particle number. At small q 2 /oo, the behavior of the dynamic 
structure factor is 



S(q,u) 



(27) 



457T 2 LO 7 / 2 

The dynamic structure function has a singularity near u = q 2 /2, which is the single-particle 
peak. It would be incorrect to take the limit x = q 2 /2u — > 1 in our formula at fixed q or fixed 
u>. Our result is strictly valid in the regime of u — > oo, q — > oo, fixed x = q 2 /2u ^ 1. The 
weaker singularity at u = q 2 /4 (x = 2) is associated with the fact that near this point, the 
two particles that are kicked out have small relative momentum, and thus interact strongly 
with each other. 



V. RESPONSE OF UNITARY GAS TO RAPID CHANGES OF THE INVERSE 

SCATTERING LENGTH 



We consider in this section the following problem. Assume we make the inverse scattering 
length a = a -1 exhibit periodic behavior in time, 



a{t) = a (t) = «o cos tut, 



(28) 



and assume the amplitude «o to be small, and the frequency of the oscillations to be large 
compared to the Fermi energy, u ^> ep. The question is: at what rate is energy deposited 
into the system? 

We know that the inverse scattering length is coupled to the operator 0*0, i.e., introducing 
a finite inverse scattering length corresponds to adding into the Lagrangian a term 



(29) 
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The energy deposited into the system, in unit volume and in unit time, can be computed 
from the formula 

d i = Sr^ lmGMl (30) 

where Goo is the Greens function of the the operator O = 0*0. The OPE of two O operators 
is easy to compute, 

-i{T 0*0(t, x) 0*0(0)) = (0*0) [D(t, x) + D(-t, -x)] + • • • , (31) 

where £)(£, x) is the scalar propagator [Eq. (5)]. From this we find ImGoo(w,0) = 
— 47rw~ 1 / 2 (0*0). Thus we find the rate of energy deposition, 

i = s<™^- (32) 

If a(t) has the form of a sudden pulse, i.e., is nonzero only in a finite time interval, 
the total energy deposited into the system, per unit volume, will be related to the Fourier 
transform a(u>) = f dte lLx)t a(t) by 

(0*0) [°°duj 



2tt 



^v^laMI 2 . (33) 



VI. CONCLUSION 



In this paper we have shown that the structure of a unitary Fermi gas, at short distance 
and length scales, can be obtained from the operator product expansion. We computed the 
OPE of the density operator at two points and derived an expression for the dynamic struc- 
ture factor. The same technique is applied to the problem of finding the energy deposited 
into a unitary Fermi gas by a rapid, small oscillation of the scattering length. 

As mentioned above, the two-body contribution to S'(q, cu) vanishes for u < q 2 /A: this is 
the two-body threshold. If u> > q 2 /6, it is possible to transfer the energy and momentum to 
three particles instead. In general, when n— 1 < x < n, where x is the "Bjorken x" variable, 
the dominant contribution to the dynamic structure factor S(u, q) in the large u, fixed x 
regime is due to the lowest dimensional operator in the OPE that contains n creation and 
n annihilation operators. The dimension of this operator is 2A n , where A n is the ground 
state energy of a system of n particles in a harmonic potential. Thus we find 

S(q,u;)~<OtO n )J^, (34) 

in the limit u — > oo, x = fixed, and n — 1 < x < n (n > 3). For example, the n = 3 
body contribution to the dynamic structure factor is down by (j~ 3 - 77272 since A 3 ps 4.27272. 
The larger n is, the faster the contribution of n-body physics decreases with increasing 
momentum. 
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The operator product expansion cleanly separates "hard" (large-momentum and/or large- 
energy) physics from soft physics. The computation of the OPE coefficients involves com- 
puting diagrams in vacuum, with a finite number of particles in the intermediate state. Thus 
the OPE is a natural way of separating few-body from many-body physics. It application 
to cold atom physics should be further investigated. 
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